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_^ , Abstract. This is a review of the main ideas of the inverse scattering method 

(ISM) for solving nonlinear evolution equations (NLEE), known as soli- 
ton equations. As a basic tool we use the fundamental analytic solutions 

'rys • x^ {^7 ^) of the Lax operator L{X). Then the inverse scattering problem for 

L(A) reduces to a Riemann-Hilbert problem. Such construction has been ap- 
plied to wide class of Lax operators, related to the simple Lie algebras. We 
construct the kernel of the resolvent of L{X) in terms of x^ {x, A) and derive 
the spectral decompositions of L(A). Thus we can solve the relevant classes 
of NLEE which include the NLS eq. and its multi-component generaliza- 
tions, the A^-wave equations etc. Applying the dressing method of Zakharov 

tH- I and Shabat we derive the A^-soliton solutions of these equations. 

^D ' Next we explain that the ISM is a natural generalization of the Fourier trans- 

form method. As appropriate generalizations of the usual exponential func- 
tion we use the so-called "squared solutions" which are constructed again 

^^ , in terms of x^{x, A) and the Cartan-Weyl basis of the relevant Lie algebra. 

f^ ' One can prove the completeness relations for the "squared solutions" which 

in fact provide the spectral decompositions of the recursion operator A. 
These decompositions can be used to derive all fundamental properties of the 
^ . corresponding NLEE in terms of A: i) the explicit form of the class of in- 

^ I tegrable NLEE; ii) the generating functionals of integrals of motion; iii) the 

• ^ . hierarchies of Hamiltonian structures. We outline the importance of the clas- 

rN ' sical i?-matrices for extracting the involutive integrals of motion. 



1. Introduction 

The modern development of the soliton theory in the last three decades of the 20- 
th century has lead to a number of important applications and developments in 
several areas of contemporary physics and mathematics, see [41, 12, 6, 3, 2]. In 
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this review I will outline the basic ideas of the inverse scattering method (ISM) on 

the example of the nonlinear Schrodinger equation (NLS) and its multicomponent 

generaUzations. 

The integrability of the well known (scalar) NLS eq.: 

iqt + qxx + 2e\q{x,t)fq{x,t)=0, e = ±1. (1) 

was discovered by Zakharov and Shabat in their pioneer paper [45] which strongly 
stimulated the search of other important integrable nonlinear evolution equations 
(NLEE). After the Korteweg-de Vries equation, this was the second NLEE inte- 
grable by the ISM. In the next few years the number of new integrable NLEE was 
growing quickly: the modified KdV eq. [39], the A^-wave equations [40, 29, 31], 
the vector NLS [34], the Toda chain [35], the principal chiral field eq. [43, 42] etc. 

The simplest nontrivial multicomponent generalizations of NLS is the vector NLS 
eq. known also as the Manakov model [34] : 

iqt + qxx + 2{q^q)q{x, t) = 0, (2) 

where g is an n-component complex-valued vector: 

/ qi{x,t) 
q{x, t) = : 

\Qnix,t) 

tending to zero fast enough for x -^ itcx). Another version of the Manakov model 
is: 

m + Qxx + 2{q^eq)q{x, t) = 0, (3) 

where e = diag (ei, . . . , e„) and e„ = ±1. Both the scalar and vector NLS equa- 
tions (especially the one with n = 2) find wide applications in nonlinear optics, 
plasma physics etc. 

Equations (2) and (3) are particular cases of the matrix NLS eq. which is obtained 
from the system: 

^Qt + (lxx + '^Qrq{x,t) =0, 

-ir-t + r^o; + 2rqr{x, t) = 0. (4) 

with appropriate reductions (involution). Here q{x,t) and r'^{x^t) are n x m- 
matrix- valued functions of x and t with n > 1, m > 1 which are smooth enough 
and tend to zero fast enough for x -^ ±oo. The best known involution compatible 
with the evolution of (4) is 

r = B_q^B^\ i?± = diag (7?f , . . . ,7?±), {r]tf = l. (5) 

and the corresponding MNLS equation is of the form: 

iQt + Qxx + 2qB^q^B^^q{x, t) = (6) 
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For n = m = 1 and r = eq* the system (4) goes into the scalar NLS equation; for 
m = 1 and n > 1 and with appropriate choice of the involution (5) eq. (4) can be 
transferred into the Manakov model or into eq. (6). 

The MNLS (6) is known to be closely related to the symmetric spaces [13]. All 
these versions of NLS are solvable by applying the ISM. Their Lax representation 
[L, M] = is provided by a generalization of the Zakharov-Shabat system: 

LtP = A— + Qix, t) - xA i;{x, A) = 0, (7) 

M^ = (i^ + yo(x, t) + A^i (x, t) - 2X^A i^ix, A) = 0, (8) 

where Q{x,t) and J are (n + tti) x (n + m) matrices with compatible block 
structure and Vo(x, t), Vi{x, t) are expressed in terms of Q and its x-derivative: 

Vi{x,t) = 2Q{x,t), Vo{x,t) = -[Q,adfQ] + 2mdfQ^. (10) 

Remark 1. One can consider Lax operators which are more general than (7). 
For example, like in [40, 11, 15] one can choose Q{x, t) and J to be elements of 
a simple Lie algebra q such that J G f) and Q{x,t) = [J,Q{x,t)]. Such form 
of Qix, t) can always be achieved by a gauge transformation commuting with J. 
Such Q{x, t) span the co-adjoint orbit of q passing through J and can be viewed 
as the tangent plane to the homogeneous space Q j J . Here Q is the Lie group with 
Lie algebra g and J is the subgroup ofQ commuting with J. Our choice ofQ{x, t) 
in eq. (7) corresponds to the symmetric space SU{n + m)/S{U{n) (g) U{ra)), see 
also [13]. 

Remark 2. An effective tool to impose involutions is the reduction group intro- 
duced by Mikhailov [36]. The involution (5) (or 7j2-reduction) can be written as: 

BU\x, t, X*)B-^ = U{x, t, A), (11) 

where B is an automorphism ofg matrix such that B^ = 1, [J, S] = and: 

U{x,t,X)=Q{x,t)-XJ, B=(^^+^_y (12) 

Reductions leading to new types of MNLS systems are demonstrated in [21, 22]. 

Section 2 is devoted to the direct scattering problem for the Lax operator L. Our 
analysis is based on the notion of fundamental analytic solution (FAS) which al- 
lows one to prove that the scattering problem for L is equivalent to a Riemann- 
Hilbert problem (RHP) for FAS. Two minimal sets of scattering data 7i 2 are in- 
troduced and shown to determine uniquely both the scattering matrix T(A) and the 
corresponding potential Q{x). 
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In Section 3 we approach the solution of the inverse scattering problem (ISP) 
through the RHP. The dressing Zakharov-Shabat method [44, 43] for the sym- 
metric spaces SU{n + m)/S{U{n) ® U{m)) is outlined and used to construct 
explicitly singular solutions of the RHP. Thus we derive reflectionless potentials 
for L and soliton solutions for the relevant NLEE. We also define the resolvent of 
L through the FAS and prove the completeness relation for the Jost solutions of L. 

Section 4 starts with the Wronskian relations as a tool to study the mapping be- 
tween the potential Q{x) and the scattering data of L generalizing the results of 
[7, 8, 1, 23]. Using them one is able to introduce the sets of 'squared solutions' 
{*} and {$} and prove that they are complete in the space of all allowed po- 
tentials. This makes more precise and explicit the results of [18]. We derive the 
expansions of Q{x) and its variation ad J^(5(5(x) over {*} and {$} and demon- 
strate that the elements of the minimal sets of scattering data 7i 2 and their varia- 
tions appear as expansion coefficients. We also introduce the generating operators 
A-t for which {*} and {$} are sets of eigen- and adjoint functions. 
The tools developed in Section 4 are used in Section 5 to describe the fundamental 
properties of the NLEE. There we prove a theorem stating the equivalence of the 
NLEE to a corresponding set of linear evolution equations for the scattering data. 
Next we derive the hierarchy of the integrals of motion from the principle series in 
terms of Q and show that they play the role of Hamiltonians for the MNLS type 
equations. We display the hierarchy of Hamiltonian structures for these NLEE. In 
short we have demonstrated the complete analogy between the ISM and the usual 
Fourier transform thus generalizing the results of [1, 30, 32, 15, 26, 27]. 

Using the method of the classical i?-matrix [12] we derive the Poisson brackets 
between the elements of the scattering matrix. As a consequence we prove that the 
integrals of motion from the principle series are in involution. 

We end by a brief discussion in Section 6 on related methods and topics. 

2. Direct and inverse scattering problems for L 
2.1. The scattering problem for L 

Here we briefly outline the scattering problem for the system (7) for the class of 

potentials Q{x, t) satisfying the following: 

Condition CI: Q{x,t) are smooth enough and fall off to zero fast enough for 

X -^ ±00 for all t. 

Condition C2: Q{x,t) is such that L has a finite number of simple eigenvalues 

A=^ G C± for all t. 

In this subsection t plays the role of an additional parameter; for the sake of brevity 

the t-dependence is not always shown. Condition C2 can not be formulated as a 
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set of explicit conditions on Q{x, t); its precise meaning will become clear below. 
The main tool here are the Jost solutions defined by their asymptotics at x -^ ±00: 

lim V(x,A)e'^^^ = II, lim (/.(x, A)e^^-^^ = II, (13) 

Along with the Jost solutions we introduce 

i{x, A) = V(x, A)e^^•^^ ^{x, A) = 0(x, A)e'^^"; (14) 

which satisfy the following linear integral equations 

e(x,A) = II +i rdye-i^-^(--^)Q(y)e(y,A)e'"^(^-^), (15) 

J 00 

(^(x, A) = II + i r dye"'^^(--^)Q(y)<^(y, A)e'"^(^-^) . (16) 

J —00 

These are Volterra type equations which, as is well known always have solutions 
providing one can ensure the convergence of the integrals in the right hand side. 
For A real the exponential factors in (15) and (16) are just oscillating and the con- 
vergence is ensured by condition CI. 

Obviously the Jost solutions as whole can not be extended for Im A 7^ 0. However 
some of their columns can be extended for A G C+, others - for A € C_. Indeed, 
the equation (15) for the first column of ^(x. A) contains only the exponential factor 
giA(a:-j/) yyhich falls off for Im A < 0. More precisely we can write down the Jost 
solutions ^(x. A) and (f){x, A) in the following block-matrix form: 

V^(x,A) = (!V-(x,A)),|V+(x,A))) , </.(x,A) = (!</>+ (x. A)), \r{x,\))) , 

(17) 
where the superscript + and (resp. — ) shows that the corresponding block-matrix 
allows analytic extension for A G C+ (resp. A € C_). 

Solving the direct scattering problem means given the potential Q{x) to find the 
scattering matrix T(A). By definition T(A) relates the two Jost solutions: 

,^(x,A) = ^(x,A)r(A), T(A) = (^+[^1 ~a-{\)) ^^^^ 

and has compatible block-matrix structure. In what follows we will need also the 
inverse of the scattering matrix: 

^{x, A) = 0(x, A)f(A), f(A) ^ ( _^^i^]^ ^;|^J ) , (19) 

where 

c-(A) = a+(A)(ll + p-p+)-^ = (II + r+T-)-^a+(A), (20a) 

d-{\) = a+{X)p'{X){t + p+p-)-^ = (II + r+T-)"V+(A)a-(A), (20b) 
c+(A) = d-(A)(ll + p^p'Y^ = (II + T-T+)-^a-{\), (20c) 
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d+(A) = d-(A)p+(A)(ll + p-p+)-^ = (11 + r-T+)-V-(A)a+(A). (20d) 

The diagonal blocks of both T(A) and T(A) allow analytic continuation off the real 
axis, namely a+(A), c+(A) are analytic functions of A for A € C±, while a^(A), 
c^(A) are analytic functions of A for A € C±. 

By p^(A) and r^(A) above we have denoted the multicomponent generalizations 
of the reflection coefficients (for the scalar case, see [1, 7, 23, 32]): 

p^{X) = b^a^{X) = c±d±(A), r±(A) = a±fo^(A) = d^c^{X), (21) 

We will need also the asymptotics for A ^ oo: 

lim 0(x,A)e'^^^ = lim V(x,A)e'^^^ = II, lim r(A) = II, (22) 

A^— oo A— too A^oo 

lim a+(A) = lim c~(A) = 11, lim a" (A) = lim c+(A) = II. 

A^oo A— >oo A— >oo A-^oo 

The inverse to the Jost solutions ip{x, A) and (j){x, A) are solutions to: 

i^-V'(a;,A)(Q(x)-AJ)=0, (23) 

satisfying the conditions: 

lim e"'^-^''V'(2^,A) = II, lim e"'^'^''<^(x. A) = II. (24) 



Now it is the collections of rows of ^(x. A) and (j){x, X) that possess analytic prop- 
erties in A: 

Just like the Jost solutions, their inverse (25) are solutions to linear equations (23) 
with regular boundary conditions (24); therefore they can have no singularities 
in their regions of analyticity. The same holds true also for the scattering matrix 

T(A) = ip{x, \)(t){x, A) and its inverse T(A) = (/)(x, \)ip{x^ A), i.e. 

a+(A) = (V'+(x,A)|,/.+ (x,A)), a-(A) = (V'-(x,A)|r(x,A)), (26) 

as well as 

c+(A) = (<^+(x,A)|^+(x,A)), c-(A) = (,^-(x,A)|V~(x,A)), (27) 

are analytic for A € C± and have no singularities in their regions of analyticity. 
However they may become degenerate (i.e., their determinants may vanish) for 
some values X^ G C± of A. Below we analyze the structure of these degeneracies. 
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2.2. The fundamental analytic solutions 

The next step is to construct the fundamental analytic solutions of (7). In our case 
this is done simply by combining the blocks of Jost solutions with the same analytic 
properties: 

X+(x,A) ^ {\cP+),\^^)) {x,X) = cPix,X)S+{X) = ^ix,X)T-{X), 

X"(x, A) = (IV^-), ir )) (x, A) = ct>{x, X)S-{X) = i;{x, A)T+(A), (28) 

where the block-triangular functions S^{X) and T^(A) are given by: 

These triangular factors can be viewed also as generalized Gauss decompositions 
(see [28]) of T(A) and its inverse: 

r(A) = T-(A)5^(A) =r+(A)5"(A), 

f(A) = 5+(A)T"(A) = 5-(A)t^(A). (30) 

The relations between c^(A), d (X) and a^(A), & (A) in eq. (20) ensure that 
equations (30) become identities. From eqs. (28), (29) we derive: 

X+(x,A) = x"(x,A)Go(A), Go(A) = ^-(A)(ll + K-(A)), (31) 

X-(x,A) = x+(x,A)G'o(A), Go{X) = D+{X){i-K+{X)), (32) 

valid for A € M, where 

Obviously the block-diagonal factors -D^(A) and D^{X) are matrix- valued ana- 
lytic functions for A G C±. Another well known fact about the FAS x^(x, A) 
concerns their asymptotic behavior for A ^ ±oo, namely: 

lim X^{x,X) = 11, (35) 

A— >oo 

where we have introduced: 

X±(x,A) = x^(x,A)e'^^^ (36) 

In the derivations that follow the analyticity properties of X^{x, A) for A € C± 
and eq. (35) will play crucial role. 
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2.3. The Gel'fand-Levitan-Marchenko equation 

The first method to solve the ISM is based on the Gerfand-Levitan-Marchenko 
(GLM) equation. For the Zakharov-Shabat system it is well known; for the block- 
matrix case - see the pedagogical exposition in [2]. The GLM equation is an in- 
tegral equation for the transformation operator /C which relate the Jost solution 

ip{x, A) and the 'plane waves' exp(— iA Jx) as follows: 

/•oo 

V'(x,A)e^^^^ = II + / ds/C(x,s)e^^-^(^-*\ (37) 

Jx 

Then the GLM equation is a Volterra-type integral equation of the form: 

/C(rc, y) =T{x + y)+ ds /C(x, s)T{s + y), (38) 

Jx 

where the kernel J-'{x) is expressed in terms of the scattering data of the operator 
L as follows: 

/9^(A) exp(— iAx) 

/9+(A)exp(iAa;) 

pj exp{-iXjx) \ 

p+ exp(iA+x) J ' ^^^^ 

where p^(A) are the reflection coefficients introduced in (21) above, X^ are the 
discrete eigenvalues of L and the constants p^ characterize the norming of the Jost 
solutions ^^(x, A^), see [2]. 

Let us now explain how the GLM eq. allows one to solve the ISP for L. Indeed, 
given the scattering matrix T(A) we can construct the reflection coefficients /5^(A), 
the coefficients p^ and the discrete eigenvalues A^, j = 1, . . . , A^. This provides 
us with the kernel J^{x) of the GLM eq. Since this is a Volterra type equation it 
always has a solution. If one is able to solve it and construct the transformation 
operator /C(x, y) then the corresponding potential Q{x) can be recovered through: 

Q{x) = iJ[J,}C{x,x)]. (40) 

Thus, given the scattering data we recover the corresponding potential Q{x). 

It is well known that in the reflectionless case p^ = p" = the kernel J^{x) be- 
comes degenerated and the GLM equation can be solved exactly thus providing the 
reflectionless potentials for L, which in turn are directly related to the soliton so- 
lutions of the corresponding NLEE. We will derive them below using the method, 
known now as the dressing method [34, 44, 42]. 
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2.4, Reductions of L. 

Typically one goes from the system (4) to the MNLS equations (1) with e = ±1 
by imposing the condition r = eq^ which is known as reduction condition. The 
theory for constructing such reductions was proposed by A. Mikhailov in [36] and 
developed further in [19, 20, 33]. 
Here we analyze just Z2-reduction, which are most widely used in the literature: 

B{U\x,t,\*))B-^ = U{x,t,\), U{x,t,X) = Q{x,t)-XJ, B^ = t. 

(41) 
Choosing B to be constant block-diagonal matrix: 

we find that eq. (41) leads to: 

r = B^q^B^^, q = B+r^BZ^, (43) 

More precisely the reduction (41) means that: 

{x+{x,t,X*))^ = B-\x'~{x,t,X))-'B, (44) 

where x^ and x'' are conveniently chosen solutions of eq. (7). If we identify x^ 
as the FAS (28) analytic for A G C± then 

(x+(x,t,A*))t = B-\x-{x,t,X)D-{X))-'B. (45) 

This result is derived easily by comparing the asymptotics of both sides of eq. (45) 
for X -^ iboo. Skipping the detail we will list here the basic consequences of the 
reductions for the scattering data of L. For the scattering matrix we get: 

T^{t,X*) = Bf{t,X)B-^, (46) 

or in 'components' (skipping the t-dependence): 

(a±(A*))t = B±c^{X)B±, (b^(A*))t = B±d^{X)B^, A G M, 

(p-(A))t = B^p+{X)B+, (T-(A))t = i?+T+(A)i?_, (47) 

{D+{X*))^ = BD-{X)B. 

For other interesting choices leading to new reductions of A^-wave and MNLS -type 
equations, see [19, 20, 21, 22]. 

2.5, The Riemann-Hilbert problem 

The eqs. (31) and (32) can be written down as: 

X+{x,X)=X-{x,X)D-{X){i + K-{x,X)), XgR (48) 
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Figure 1. The contours j± = R U 7±c 



or 



X-{x,X) =X+(x,A)i)+(A)(]l -K+{x,X)), 



where 



K^{x,X) 



-iXJx 



K^iX)e 



iXJx 



AG 



(49) 



(50) 



Eq. (48) (resp. eq. (49)) combined with (35) is known in the Uterature [14] 
as a Riemann-Hilbert problem (RHP) with canonical normalization. It is well 
known that RHP with canonical normalization has unique regular solution; the 
matrix-valued solutions Xq{x, A) and Xq{x, A) of (48), (35) is called regular if 
det Xq (x, A) does not vanish for any A G C±. 

Let us now apply the contour-integration method to derive the integral decomposi- 
tions of X^{x, A). To this end we consider the contour integral: 



^i(A) 



1 



d^ 



■X+{x,ii) 



1 



d^ 



27ri J^^ fi — X 27ri /■y_ fi — X 

where A G C+ and the contours 7-1- are shown on fig. 1. 
First we evaluate c7i(A) by Cauchy theorem. The result is: 



X-{x,fi)D-ifi), (51) 



N 



Ji(A) =X+{x,X)+J2 ^^^ 



x-{x,^i)D-{^i) 

fi — X 



(52) 



We can also evaluate Ji(A) by integrating along the contours. In integrating along 
the infinite semi-circles of 7±,oo we use the asymptotic behavior of X^{x, A) for 
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A ^ oo. The result is: 



1 r°° du 

MX) = 11 + — / -J^x-{x,fi)D-{fi)K-{fi), (53) 

where in evaluating the integrand we made use of eq. (48). Equating the right hand 
sides of (52) and (53) we get the following integral decomposition for X+(x, A): 

X+(.,A) = Il-f Res ""'^^'^^^'^^^ 

1 1"°° du 
+ — / -J^X-{x,fi)D-{fi)K-{fi), (54) 

ivri j-oo /U — A 

Quite analogously we derive the decomposition for X^{x, A): 

X-(.,A) = a-ERes^"(-'^)f(^) 

+ ^ r ^^+(^,/^)^+(a*)^+(a*), (55) 

zin J-oo /U — A 

Equations (54), (55) can be viewed as a set of singular integral equations which are 
equivalent to the RHR For the MNLS these were first derived in [34]. 



2.6. The minimal set of scattering data 

Obviously, given the potential Q{x) one can solve the integral equations for the 
Jost solutions which determine them uniquely. The Jost solutions in turn determine 
uniquely the scattering matrix T(A) and its inverse T(A). But Q{x) contains at 
most 2mn independent complex- valued functions of x. Thus it is natural to expect 
that at most 2mn of the coefficients in T(A) for A € M will be independent; the 
rest must be functions of those. 

The set of independent coefficients of T(A) are known as the minimal set of scat- 
tering data. As such we may use any of the following two sets % = 7^ c U 7^ d^ 

Ti,c^{/9+(A),p-(A), A€M}, ^i,d^{pf,Af}J^^, 

T2,c^{r+(A),r-(A), A G M} , ^M = {^f , A±}J^^ , (56) 

where the reflection coefficients /5^(A) and r^(A) were introduced in eq. (20), 
A^ are (simple) discrete eigenvalues of L and p^ and r^ characterize the norming 
constants of the corresponding Jost solutions. 
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The reflection coefficients p^(A) and r^(A) are defined only on the real A-axis, 
while the diagonal blocks a^(A) and c^(A) (or, equivalently, D^{X)) allow ana- 
lytic extensions for A € C-t-. From the equations (20) there follows that: 

a+(A)c-(A) = (II + p-p+{X))-\ a-(A)c+(A) = (II + p+p-(A))-\(57) 

c-(A)a+(A) = (II + T+T-{X))-\ c+(A)a-(A) = (II + r-T+(A))'^(58) 

Given 7i (resp., T2) we determine the right hand sides of (57) (resp. (58)) for 
A G M. Combined with the facts about the limits: 

lim a+(A) = II, lim c"(A) = II, 

A^oo A^oo 

lim a" (A) = 11, lim c+(A) = II, (59) 

A^oo A— >oo 

each of the relations (57), (58) can be viewed as a RHP with canonical normaliza- 
tion. Such RHP can be solved explicitly in the one-component case (provided we 
know the locations of their zeroes) by using the Plemelj-Sokhotsky formulae [14]. 
These zeroes are in fact the discrete eigenvalues of L. One possibility to make use 
of these facts is to take log of the determinants of both sides of (57) getting: 

A+(A) + C-(A) = -lndet(ll + p-p+(A), A G M, (60) 

where 

A^{X) = Indet a^(A), C^(A) = Indet c^(A). (61) 

Then Plemelj-Sokhotsky formulae allows us to recover ^^(A) and C^(A): 

-00 p 



\ r°° c\n ^ A - A+ 

2tt J -00 /U - A ^ A - Aj- 

where ^(A) = ^+(A) for A G C+ and ^(A) = -C-(A) for A G C_. In deriving 
(62) we have also assumed that A are simple zeroes of ^^(A) and C^(A). 

Remark 3. The sets (56) were first derived for m = n = 1 in [30], see also [23]. 
Here we are using the gauge covariant formulation developed in [26]. 

If we impose the reduction condition (47) then we get the following constraints on 
the sets Ti 2: 



p-{X) = (S_p+(A)i?+)t, 


P, 


= (i?-P+5+)t, 


V 


= (A/)*, 


(63) 


T-{X) = {B+T+{X)B^)\ 


^i 


= (B+r;i?_)t, 


^^ 


= (\^)*' 


(64) 


where j = 1, . . . , A^, 













Remark 4. For certain choices of the reduction conditions (40), such as Q = 
— Q' the generalized Zakharov-Shabat system L{X)'tp = can be written down as 
an eigenvalue problem Cip = Xip{x, A) where C is a self-adjoint operator. The 
continuous spectrum of C fills up the whole real X-axis thus 'leaving no space ' 
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for discrete eigenvalues. Such Lax operators have no discrete spectrum and the 
corresponding MNLS do not have soliton solutions. 

The full RHP's (57), (58) for n and tti > 1 do not allow explicit solutions. How- 
ever, from the general theory of RHP [14] one may conclude that (57), (58) allow 
unique solutions provided the number and types of the zeroes A^ are properly cho- 
sen. 

Thus we can outline a procedure which allows one to reconstruct not only T(A) 
and T(A) and the corresponding potential Q{x) from each of the sets Ti,i = l, 2: 

i) Given Ti (resp. 7^) solve the RHP (57) (resp. (58)) and construct a^(A) 

andc±(A)for A G C±. 
ii) Given 7i we determine b (A) and d {X) as: 

6±(A) = />±(A)a±(A), d±(A) = c±(A)/>±(A), (65) 

or if 72 is known then 

b±(A) = a±(A)r±(A), d^{X) = r±(A)c±(A). (66) 

iii) The potential Q{x) can be recovered from 7i by solving the GLM eq. and 
using eq. (39). 

Another method for reconstructing Q{x) from Tj uses the interpretation of the 
ISM as generalized Fourier transform, see [1, 32, 30, 23, 15, 27, 26] and Section 4 
below. 



3. The RHP and the inverse scattering problem 

There are two pioneer results that enhanced the development of the ISM. The first 
is the existence and explicit construction of the FAS for rather general class of Lax 
operators [37, 38, 40, 4, 5, 27]. The second is the discovery of the equivalence of 
the ISP for L to a RHP for the FAS [38, 44, 43, 42]. 

Here we are considering rather special Lax operators for which the construction 
of the FAS does not require special efforts. Here we first prove the equivalence 
between RHP and the spectral problem L. Next we show how the FAS can be used 
to construct the resolvent of L and for the analysis of its spectral properties. 

The use of the RHP rather than the GLM equations allowed for important effec- 
tive extensions of the ISM; it allowed the treatment of generic Lax operators with 
rational dependence on A [43, 42, 33], as well as of ones related to simple Lie al- 
gebras [15, 19, 20, 21, 16, 27]. In the next subsections we will make use of these 
developments to describe two types of soliton solutions for the MNLS. 
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3.1. Equivalence of RHP to ISP 

The next important step is the possibility to reduce the solution of the ISP for the 
generalized Zakharov-Shabat system to a (local) RHP. Indeed the relation (31) can 
be rewritten as: 

X+{x,t,X) = X-{x,t,X)Gix,t,X), AgM, (67a) 

G{x,t,X) =e-'(^-^^-^W*)Go(A)ei(^^^--^(^)*\ (67b) 

Go(A) = 5-(A,t)5+(A,t)[^^; (67c) 

in other words the sewing function G{x, t, A) satisfies the equations: 

fir* r\C 

i—-X[J,Gix,t,X)]=0, i— + [f(X),Gix,t,X)]=0, (68) 

Here /(A) € f) determines the dispersion law of the NLEE. 

Theorem 1 ([38]). Let X~^{x, t, A) and X~{x, t, A) be solutions to the RHP (67) 
with canonical normalization (35) allowing analytic extension in Xfor X € C± 
respectively. Then x^{x^ t, A) = X^{x, t, X)e^'^'^^ are fundamental analytic solu- 
tions of both operators L and M, i.e. satisfy eqs. (7), (8) with 

Q{x, t) = lim A ( J - X^{x, t, X).JX^{x, t, X)) . (69) 

Proof: Let us assume that X^(x, t, A) are regular solutions to the RHP and let us 
introduce the function: 

dX-^ 
g^(x,t,X) = i——X^(x,t,X) + XX'^(x,t,X)JX^(x,t,X). (70) 
ax 

If X^{x, t, A) are regular then neither X^{x, t, A) nor their inverse X^{x, t, A) 
have singularities in their regions of analyticity A G C±. Then the functions 
g^{x, t, A) also will be regular for all A € C±. Besides: 

lim g^{x,t,X) = lim g^(x,t,X) = XJ. (71) 

A— +00 A— +00 

The crucial step in the proof [44] is based on the chain of relations: 

g+(x, t, A) ^=^ ili^-^GX- {x, t, A) + XX-GJGX- {x, t, A) 
ax 
dX^ / Ar~< \ 

= i^—X- (x, t, A) + X- i^G + AG JG X'ix, t, A) 
dx \ dx J 

(68) . dX ^_ ^^ ^ ^_ / ^^ ^ AG JG) X~ (x, t, A) 
dx \ / 

= i^—X~ (x,t, X) + XX- JX- (x,t, X) 
dx 

= g-{x,t,X), AgM. (72) 
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Thus we conclude that g^{x,t,X) = g~{x,t,X) is a function analytic in the 
whole complex A-plane except in the neighborhood of A ^ oo where g^{x,t,X) 
tends to A J, (71). Next from Liouville theorem we conclude that the difference 
g~^{x, t, A) — AJ is a constant with respect to A; if we denote this 'constant' by 
—Q{x,t) we get: 

g^{x,t,X)- XJ = -Q{x,t). (73) 

Using the definition oi g^{x, t, A) (70) we find that X^{x, t, A) satisfy 

i^^ + Q{x, t)X^{x, t, A) - X[J, X^{x, t, A)] = 0, (74) 

i.e. that x^{x, t, A) is a fundamental solution to L. The relation between Q{x, t) 
and X^{x, t, A) (69) is obtained by taking the limit of the left-hand sides of (73) 
for A ^ oo and using the asymptotic expansions of X^{x, A) and X^{x, A): 

oo 

X^ix,t,X) = i + J2X-'Xt{x,t), 

s=l 

oo 

X^{x,t,X) = ][ + Y,^''^tix,t), (75) 

s=l 

we get: 

Q{x,t) = lim A(J-X±jl±(x,t,A)] = [J,Xf (x,t)]. (76) 

A— »oo 

Here we also used the fact that Xi{x, t) = —X^{x, t), see eq. (75). 
Arguments along the same line applied to the functions /i^(x, t, A) 

h^{x, t, A) = i——X^{x, t, X) + 2X^X^{x, t, X)JX^{x, t, A), (77) 

can be used to prove that x^(x,t, A) are fundamental solutions also of the op- 
erator M (8). Indeed, repeating the above arguments for h^{x,t,X) we find 
that /i+(x,t, A) = h~{x,t,X) is a function analytic everywhere in C except at 
A ^ cx) where it behaves like 2A^J. From Liouville theorem it follows that 
h^{x, t, A) — 2X^J is a linear function of A equal to —Vo{x, t) — AVi(x, t). Thus: 

AY^ 

\^^ + {Vo{x, t) + XVi{x, t))X^{x, t, A) - [2AV, X^{x, t, A)] = 0, (78) 

To conclude the proof we have to account for possible zeroes and pole singularities 
of X^(x, t, A) at the points A^. Below we derive the structure of these singularities 
and show that they do not influence the functions g^{x,t, X) and h^{x,t, X). □ 



16 Vladimir S. Gerdjil<ov 



3.2. The dressing Zakharov-Shabat method for the symmetric spaces 

SU{n + m)/S{U{n) U{m)). 

Let us outline how one can, starting from a given regular solutions Xq (x,t, A) 
of the RHP, construct new singular solutions X^{x, t, A) having zeroes (singular- 
ities) at the prescribed points A^ G C±. The structure of these singularities are 
determined by the dressing factor Uj{x, t, A): 

^^{x,t,X)=Uj{x,t,X)^^{x,t,X)w-X{X), (79) 

which in our case has a simple fraction-linear dependence on A: 



A-A+ 
Uj{x,t,X) = 11 + (c,(A) - l)Pj{x,t), Cj{X) = j^- (80) 

Here Pj{x,t) is a projector P^ = Pj, Wj^+{X) = 11 and Wj^^{X) will be defined 

below. 

It is well known that, if Q{x, t) — XJ takes values in the Lie algebra g then 

X^{x,t,X) must take values in the corresponding Lie group ©. Therefore the 

dressing factor Uj{x, t, A) must also be element of the same group. But assuming 

Pj{x, t) is a projector we can derive that 

det u{x,t,X) = {cj{X)Y\ rj = mnkPj{x,t). (81) 

Obviously the ansatz (80) is compatible with © ~ GL{n + m), SL{n + m) 
and U{n + m), SU{n + m). The last two possibilities are realized provided 
proper reduction conditions like (63), (64) are imposed on A^ and Pj{x,t) = 

BP^{x,t)B-K 

Usually Pj{x) is chosen to be of rank 1. We will consider slightly more general 
case when rank Pj (x) = rj > I; rj is the multiplicity of the corresponding eigen- 
values A^ and rj < min(n, m). Then Pj{x) can be written in the form: 

Pj{x) = \nj) {{mj\nj))^ {mj\, (82) 

where the collections of rj bra- (resp. ket-) eigenvectors {mj\ (resp. \nj)) can be 
viewed also as rectangular {n + m) x rj (resp. rj x (n + m)) matrices. Then 
{mj\nj) will be quadratic rj x rj matrix which we assume to be nondegenerate. If 
rj = 1 eq. (82) provides the standard expression for rank 1 projector. It is easy to 
check that the left hand side of (82) satisfies identically the relation P? = Pj . 

From (79) there follows that the dressing factor u{x, t, A) satisfies the equation: 

i- — h Q(x, t)u(x, t, A) — u(x, t, X)Qo(x, t) — A[J, u(x, t, A)] = 0, (83) 
dx 

where QQ{x,t) (resp. Q{x,t)) is the potential related to the regular Xo{x,t,X) 
(resp. singular x^{x, t, A)) solution of the RHP. 
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Next we insert the anzats (80) and request that it holds identically with respect to A. 
To this end it is enough that eq. (80) holds true for A = A^, A ^ A~ and A ^ cx). 
The first two conditions lead to the following equations for {mj\ and \nj): 

i^ + ^W(x,t, A+)|n,) = 0, i^ + y(°)(x,t, A+)ln,) = 0, (84) 

i^ - K|t/W(x,t, AT) = 0, i^ - (m,|yW(x,t, A-) = 0,(85) 
U^'Hx,t,X)=Qo{x,t)-XJ, y(o)(x,t,A) = y(x,t,A)|Q^Q,^. (86) 

Here V^^^ {x, t, A) is obtained from V{x, t, A) (see (8), (10)) replacing Q{x, t) by 
QQ{x,t). This construction is well defined also in the case when Xo (^j-^) ^^ 
singular solutions to the RHP, provided they are regular for A = A^. In order 
to avoid technicalities in what follows we will treat only the case of one pair of 
discrete eigenvalues. From eqs. (84) there follows that: 

l^i) = xtj{x,t)\n°j), {mj\ = {m^j\xoj{x,t), Xoj{x,t) = Xo(.x,t,Xf). 

(87) 
Note that since Xo (^' ^^ solutions of the regular RHP then Xoji^y t) exist and 
are nondegenerate. The equation (69) considered for A ^ oo gives the following 
relation between Qo{x, t), Q{x, t) and P{x, t): 

Q{x,t) = Qo{x,t) + lim X{J — Uj{x,t,X)Juj{x,t,X)) 

A— »oo 

= Qo{x,t) - (A+ - Xj)[J,P,{x,t)]. (88) 

Thus starting from a given regular solution of the RHP (and related solution Qo{x, t) 
to the NLEE) we can construct a singular solution to the RHP and a new solution 
Q{x,t) of the NLEE depending on the A^ and on the eigenvectors of Pj{x). If 
we start from the trivial solution Qo{x, t) = of the NLEE then we will get the 
one-soliton solution of the NLEE. Repeating the procedure N times we can get the 
A'^-soliton solution of the NLEE. 

With the explicit formulae for Pj{x) and using (79) we can establish the relation- 
ship between the scattering data of the regular RHP and the corresponding singular 
one. The dressing factor Uj (x, t, A) is determined by the collections of constant po- 
larization vectors 

In?) = ( j^lj:| ) , Kl = (("^^,,1, K,,l) , (89) 

which can not be quite arbitrary. They should be such that the constant rj x rj 
matrices {rriQ Auq •) and {rriQ Auq •) are nondegenerate. 
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^n,j 


= Il + (c,(A)- 


-i)^A,' 


(92) 


^22,j 


= fl + (c,(A)- 


-l)^2-2,r 


(93) 



Now we can evaluate the limits of the dressing factor u{x, A) for x -^ ±oo and 
then derive the interrelations between the scattering matrices Tq{X) and Ti(A) cor- 
responding to the potentials Qo{x,t) and Qi{x,t). In what follows the elements 
of the scattering matrix Tq{X) will be denoted by the same letters as the ones of 
T(A) but with additional subscript 0. 
Using the explicit formulae (87), (28), (29) we get: 

^' = .i!Too^^(-)= ^" = (1'^' o) ' ^^(-) = (o 4) '^90) 

where obviously rankP^^ = rankP22 = ^j- Therefore the limits of the dressing 
factor are given by: 

x^-oo y U U22J J 

Next we have to determine the matrices Wj,±(A) in the right hand side of eq. (79) 
so that the asymptotics of the singular solution are compatible with (28), (29). This 
holds true if 

^,,+(A) = II, w,,.{X) ^ W{X) = H^'^ ^°^^ ] . (94) 

In addition we find the following relations between the scattering matrices: 

T(A) = u+To{X)uJ, (95) 
or 'in components': 

a+(A) = n+^^.a+(A), 6-(A) = n+_^.&o-(A)ti22,i> (96a) 

a- (A) = ao-(A)n22,,> ^^W = HW, (96b) 

c+(A) = n22^^.c+(A), d+(A) = u^2A(^^^ti,j^ (96c) 

c-(A) = c,{X)ut,^^, d-{X) = d,{X), (96d) 

The dressed solutions of the RHP are given by: 

X+(a:, A) = u{x, X)xo{x, A), x'{x, A) = u{x, X)xo{x, A)W^(A), (97) 

D+{X) = W{X)D+{X), D-iX) = Do{X)W{X), (98) 
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W{\) = 11 + (c,(A) - l)W,, W, = ( ^JJ'^- ^° j . (99) 



The explicit construction of the dressed FAS allow us to reveal the structure of the 
FAS at A ~ A^. Here we first formulate the expansions of a^(A), c^(A) and their 
inverse in the vicinities of A 

a±(A) = af + (A - Aj)d±(A) + 0((A - A^^)^), (100a) 

c±(A) = c± + (A - A±)c±(A) + 0{{\ - Xff), (100b) 



a± 



^'=(^) = ,, \±. + ^"^CA) + 0{{\ - Xf)\ (lOOd) 

(A - A ■ ) 



where 



af = (11 - PA,.)<' «f = (^^ - V)«^.^i'ij' (101^) 

c± = c+^.(ll - Pi+^.), c± = (AT - A+)Pi+ /+^., (101b) 

Let us outline the structure of the eigenspaces corresponding to the discrete eigen- 
values A^. To avoid technicalities in doing this we will assume that: i) L has no 
other discrete eigenvalues and ii) that the regular potential Q{x) is on finite sup- 
port. Due to ii) one can prove that all Jost solutions and their inverse, as well as 
T(A) and T(A) are meromorphic functions of A and can be extended to the whole 
complex A-plane. Considering equations (18), (19) at A = A^ we derive the rela- 
tions: 

\4'f{x))af = ±\^f{x))pf, pf = bfaf, (102a) 

\i^f{x))cf = ±|0±(x))rf , rf = djcf, (102b) 

where the index j means that we are taking the value of the corresponding function 
for A = A+. 

But the eigenfunctions corresponding to the discrete eigenvalues must be square 
integrable. A necessary condition for this is the requirement that these eigenfunc- 
tions have no exponentially growing terms for both a; ^ oo and x ^ — oo. These 
limits for xt{^) = X'^{x, Xf) = (|(^^(a:)), [^-""(x))) are equal to: 





lim_Xj{x) = e ^ \ ^^ ^ j , ^\hn^Xj{x) = e 

' (103) 

Since generically af ^ and cf ^ some of the columns of xt (x) will be expo- 
nentially growing and can not be interpreted as discrete eigenfunctions. However 
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from eqs. (96) it is clear that both a'j' and ct are degenerate matrices of rank n — tj 
and m — Vj respectively. Thus we find that rj linear combinations of columns of 
\4>'^{x)) and \ijj'^{x)) due to the relations 

a+|ni^^.) = 0, cj\nl^)=Q, (104) 



decrease exponentially for both x ^ oo and x ^ — oo. The eigenspace of L 
\^ is spanned by Vj 



related to A^ is spanned by r, linearly independent discrete eigenfunctions which 



can be chosen among |0^(x)nQ •) or \il)'^{x)n'^ •). These two sets of eigenfunc- 
tions satisfy linear relations which can be written compactly as: 

X+(x, A+)|no,,) ^=^11 - Pj{x))\nj{x)) = 0. (105) 

The Vj constant vectors |noj) which determine the discrete eigenspace can be 
viewed as 'polarization' vectors of the corresponding soliton solution and parame- 
trize its internal degrees of freedom. 

The eigenspace corresponding to A^ can be analyzed either independently along 
the same lines or by using the reduction properties of L, see subsection 2.4. Indeed, 
from (96b) and (96d) we get that only rj linear combinations of the columns of 
\(j)'j{x)) (resp. |^"(x))) do not have exponential growth for x ^ oo (resp. x -^ 
— oo). The corresponding 'polarization' vectors \n^ ■) are related to |noj) by: 

K,) = ^o';,>o,,), Dl- = D+{X+). (106) 

because 

The analog of eq. (105) for xj (x) is the following relation: 

XJ{x)Kj) = 0' Kj) = DU^t)\m,j). (108) 

Remark 5. Each of the eigenvalues A corresponds to rj -dimensional eigensub- 
space with 1 < rj < min(n, m). In particular for the vector NLS (the Manakov 
model) we may have only rj = 1. The reduction (41) relates the sets of polarization 
vectors by: 

B-^\ml^) = |no,,). (109) 

3.3. Reflectionless potentials and soliton solutions 

The simplest situation in which the dressing method can be applied is the one 
corresponding to vanishing potential and plane wave solution of (7): 

Qoix) = 0, x^{x,X) = e-'^J-^'^-'-'K (110) 



BASIC ASPECTS OF SOLITON THEORY 21 

The corresponding scattering matrix of course is To(A) = 11; therefore the corre- 
sponding reflection coefficients (21) are vanishing Pq{X) = 0, Tq (A) = on the 
whole real A axis. 

Applying the dressing method to (110) we obtain for the simplest reflectionless 
potential: 

Qi{x) = ^{X+~Xj)[J,Pj{x,t)]. (Ill) 

Obviously, due to (96) the corresponding reflection coefficients p^ (A) = 0, t^ (A) = 
are also vanishing on the whole real A axis. Therefore Qi{x) is the simplest 
non-trivial reflectionless potential of L. The scattering matrix for this potential is 
block-diagonal, but is not trivial: 

r.(A) = f<f' .-» ) (1.2) 

The reflectionless potentials are closely related to the soliton solutions of the cor- 
responding equation. Indeed, the derivation of the reflectionless potentials we con- 
sidered the time t as an auxiliary parameter. In order to go from the reflectionless 
potential to the solution of the corresponding NLEE all we need to do is to impose 
the correct t-dependence on the 'polarization vectors' \nQj) and (mojl- This can 
be determined from the dispersion law of the NLEE as follows: 

jnoj) ^ exp(-2i/(A+)t)|noj), {moj\ -^ (mojj exp(2i/(Aj)t). (113) 

Applying this procedure to Qi we get the one-soliton solution of the generic NLEE 
with dispersion law /(A) in the form: 

Qis(x,t) = 2(A+-A7) (114) 

/ -\n]{x,t))R~Hx,t){m]{x,t)\\ 

[\n]{x,t))R-\x,t){m]{x,t)\ j' 

Rj{x,t) = (m](x,t)|nj(x,t)) + (m|(x,t)|n^(x,t)), 

{m]{x,t)\ = (m],o|e'(^^"+'^<^^■)*^ {m]{x,t)\ = (m|_o|e"'^^^' "^'^'^^■^*\ 

where /q • = fo{Xj) and /(A) = /o(A) J is the dispersion law of the NLEE. In 
order to obtain the one-soliton solution for the MNLS we have to replace /o(A) J 
by /o,NLS = — 2A^. If we put rj = 1 and take into account the reduction (41) then 
Rj{x, t) can be written down as 

i?j(x, t) = 2Rqj cosh(2z>'ja: + 2pjt + ^o,j)) 

where 

u, = i(A7 - A+)/2, Jl, = i{f{Xj) - /(A+)), 
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iii\/2i2\ f ^T \"^i,oFo,i/ 

^iol"Oj)(Kol«0,i)> ^0,i = 2 In ,2 |„2 ^ 

^ \''''j,0l'''0,i/ 



thus reproducing the well know result for the one-soliton solution of the vector 
NLS and the rank-1 solutions of the MNLS [34, 2]. 

We can apply the dressing procedure again, starting with the potential Q\^(x^i) 
and the corresponding FAS Xfs(x, A). Doing this we have to use a dressing factor 
U2s(x, A) like in (80) but with new locations A^ for the eigenvalues and a new 
choice for the projector Pk{x) which may have different rank r^ and sets of vectors 
|nfc(x)), {mk{x)\. Choosing A^ ^ A^ the FAS xfs(^'-^) ^i^^ ^^ regular at the 
points A = A^ which is all that is required for the procedure to be valid. 

Repeating the procedure several times we get more and more complicated poten- 
tials which are still reflectionless. 

Another way to derive these potentials and the related A^-soliton solutions of the 
NLEE consists in using equations (54), (55) with K^(x, A) = K^{x, A) = 0. To 
explain this better following [34] we need to evaluate the residues of X^ {x, X)D^ (A) 
at A = A^. Using eq. (102) we obtain: 

Res X+{x,X)D+{X) = XUx)JCJ{x), (115) 

A=A+ 

Res X~{x,X)D~{X) = -Xr{x)ICj{x), (116) 

A=A- -^ -^ 

where Xf{x) = X^{x, Xf) and 

=t ±2iA±a: ^±f^\ _ ^±„T2iA±a: 



Thus the equations (54), (55) take the form: 

TV , 

X+(x,A) = t-Y.\ —Xj{x)lCj{x), (118) 

3=1 ^ ~ ^• 

^ 1 
X-(x,A) = ]l + ^- -^X+(x)/C+(x). (119) 

Note that it is legitimate to put A = A^ in eq. (118) and A = A^ in eq. (119) 
thus getting an algebraic system of equations for X (x). After solving for X (x) 
we can insert them into eqs. (54), (55) thus recovering the FAS X^{x, X) for all 
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A G C±. The corresponding reflectionless potential is given by (see eq. (76)): 

N N 

QW(x) = -J2[J,Xr(x)ICj{x)] = J2[J,X+{x)IC+{x)]. (120) 

i=i i=i 

As we shall see in the next subsection the generalized Zakharov-Shabat system 
(7) with potential Q^^\x) has 2N discrete eigenvalues located at A^ with rj- 
dimensional eigenspaces. The corresponding A^-soliton solution QNs(a;,t) is ob- 
tained from (120) by introducing the relevant t-dependence into the polarizations 
vectors \noj) and {moj\ or equivalently, in /C- (x,t) as follows: 

pf(x, t) = pf{x)e^^<^\ T^ix, t) = Tf{x)e^^<^\ (121) 

where f^- = /o(A^) and /(A) = /o(A) J is the dispersion law of the NLEE. 

3.4. Spectral properties of L 

Using the FAS one can construct the resolvent of (7): 

R^{x, y, A) = -ix^{x, A)e±(x - y)x^{y, A), (122) 

G±(z) = Amg{^e{^z)^,±e{±z)^), 

where 9{z) is the standard step-function. 

Let us consider R^{x, y, \) as the kernel of an integral operator acting on vector- 
valued functions of x as follows: 

/CO 
dyR^{x,y,X)f{y), forAGC±. (123) 

-co 

Theorem 2. Let Q{x) satisfy conditions (C.l) and (C.2) and let A be the zeroes 
ofdet a (A). Then 

1. R^{x, y, A) is an analytic function of \ for A G C± having pole singulari- 
ties at A • G C±; 

2. R {x, y, A) is a kernel of a bounded integral operator for Im A ^ 0; 

3. R {x, y, A) is uniformly bounded function for A G M and provides a kernel 
of an unbounded integral operator; 

4. R^{x, y, A) satisfy the equation: 

L{X)R^{x,y,X) = ll6{x-y). (124) 

Idea of the proof: 1. is obvious from the fact that x^i^^ A) are the FAS of 

L{X). From the definition (28) there follows that det x^ {x, A) = det a^ (A), 
i.e. x^{y, A) and consequently, R^{x, y, A) will develop pole singularities 



for all A^ for which det a^(A) = 0. 
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2. Assume that Im A > and consider the asymptotic behavior of R^{x, y, A) 
for x,y ^ oo. From equations (28), (29) we find that 

n 

R^{x,y,X) = J2X+{x,X)e-'^'(--y^e^{x-y)X+{y,X) (125) 
p=i 

Due to the fact that x^i^^ ^) has block-triangular asymptotics for x — > oo 
and A € C+ and for the correct choice of Q^{x — y) (122) we check that 
the right hand side of (125) falls off exponentially for x ^ oo and arbitrary 
choice of y. All other possibilities are treated analogously. 

3. For A G M the arguments of item 2) can not be applied because the expo- 
nentials in the right hand side of (125) Im A = only oscillate. Thus we 
conclude that R^{x, y, A) for A € M is only a bounded function and thus 
the corresponding operator R{X) is an unbounded integral operator. 

4. The proof of eq. (124) follows from the fact that L{X)x^{x, A) = and 

^^^^ = m^-y). (126) 

ax 
The theorem is proved. D 

From theorem 2, item 3) there follows that the continuous spectrum of L fills up the 
whole real A-axis with multiplicity n + m. By definition the operator L may also 
have discrete eigenvalues at the points at which R^{x, y, A) have pole singularities. 
From item 1) it follows that these are precisely the points A^. 
Let us now analyze the structure of these singularities and evaluate the correspond- 
ing residues. To this end we insert eqs. (97), (99) into (122). The result is that 
R^{x, y, A) have poles of first order in the neighborhood of A^ with residues : 

Res^ R+{x,y, A) = -i(A+ - A7)(ll - P,(x))x^,(x)e+(x - y)x+^{y)Pj{y). 
X-X. 

(127) 

Res R'{x,y, A) = i(A+ - Xj)P,{x)xo,{x)e'{x - i/)xo",(y)(ll - Pj{y)). 
x=\- 

(128) 

Formally in the right hand side of (127) there enter the discontinuous functions 

Q^{x — y). However, due to the special structure of the projectors Pj{x) (see eqs. 

(82) and (102a), (102b) we obtain the following expressions for the residues of 

i?±(x,y,A): 

Rts^R^{x,y,X)=±i\i^f{x))pf{i;f{y)\. (129) 

A . 

J 



A=A± ' '" ' ^^,^^,^,,r^,r^ 



where p^ are defined in eq. (102) and there are no discontinuities present. 
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Now we can derive the completeness relation for the eigenfunctions of the Lax 
operator (7) by applying the contour integration method (see e.g. [24, 1]) to the 
integral: 

^^^^y) = ^J dAi2+(x,y,A)-^^ dAi?-(x,y,A), (130) 

where the contours ^± are shown on the Figure 1. Skipping the details we get: 

5{x-y)J 

= ^y_^dA{|,/.+ (x,A))a+(A)(V+(y,A)| - |r (x, A))a-(A)(V-(y, A)|} 



-oo 

N 



-\y: (iV';(^))p+(V';(y)i - \^j{x))p]{^j{y)\) . 

(131) 

The completeness relation (131) is a natural generalization of the one in [24] for the 
sl{2) case. An important difference here is that now we have matrix- valued spec- 
tral functions a^(A) whose zeroes determine the location of the discrete eigenval- 
ues. 



Remark 6. When both n > 1 and m > 1 there are two possible definitions of 
simple eigenvalues. One of them used in [2] /. to define Xf as simple ifdet a±(A) 
have simple zeroes for A = A . The other possible definition is: the eigenvalues 
A are simple if the resolvent of L has simple poles at A . The singular solutions 
of the RHP (80) correspond to simple poles of the resolvent (122) although in the 
neighborhood of Xc^ A • we have that det a (A) behaves like (A — A • Y'L 

4. The Generalized Fourier Transforms for Nonregular J 

The main result in this section consists in the fact that the analysis of [15] can be 
applied also to nonregular choices of J. In addition we briefly outline also how to 
take into account the presence of discrete spectrum of L. Skipping the details of 
the proof we formulate the completeness relation for the 'squared' solutions of L. 

4.1. The Wronskian relations 

The analysis of the mapping J^: M. ^ T between the class of allowed potentials 
J\A and the scattering data of L starts with the Wronskian relations, see [7, 6] for 
s/(2)-case and [18] for the block-matrix case (7). As we shall see, they would 
allow us to 

1. To formulate the idea that the ISM is a GFT; 

2. To determine explicitly the proper generalizations of the usual exponents; 
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3. To introduce the skew-scalar product on M which provides it with a sym- 
plectic structure. 

All these ideas will be worked out for the system (7). With (7) one can associate 
the systems: 

i-^ -i){x,t,X)U{x,t,X) =0, U{x,X) = Q{x)-XJ, (132) 

i— — + 6U(x, t, X)ij(x, t, A) + U(x, t, X)6'il^(x, t,X) =0 (133) 

dx 

i-- - A J^(x, t, A) + U(x, t, X)tp(x, t, A) = (134) 

dx 

where 5tP corresponds to a given variation 6Q{x, t) of the potential, while by dot 
we denote the derivative with respect to the spectral parameter. 

We start with the identity: 

/•oo 

lOO 



ixJxix,X)-J)\Z_^ = -i / dxxmx),J]xix,X), (135) 



— CO 



where x(a;, A) can be any fundamental solution of L. For convenience we choose 
them to be the FAS introduced above. 

The l.h.side of (135) can be calculated explicitly by using the asymptotics of 
X^(x, A) for X -^ ±oo, (28), (29). It would be expressed by the matrix elements 
of the scattering matrix T(A), i.e., by the scattering data of L as follows: 

{x^Jx^ix,X) - J)|- _ = -2 (,+*J^^ %^")) (136a) 

ix~Jx-ix,X) - J)ir=-oo = -2 (/(;,) %^^^) (136b) 

We will show that these Wronskian relations allow us to express the elements of 
each of the sets %, i = 1,2 in eq. (56) as integrals from the potential Q{x) 
multiplied by some bilinear combination of eigenfunctions of L called 'squared 
solutions'. Let us now analyze the equations obtained from (135) after multiplying 
by the matrix Eat, {Eab)cd = ^ac^bd and taking the trace. Such operation will 
produce the a, 6-matrix element of eq. (135). In the right hand side of this equation 
we can use the invariance properties of the trace and rewrite it in the form: 



/oo 
dxtr([Q(x),J]e,fe(x,A))(137) 
-oo 

e 

eab{x, X) = xEabX{x, A), e„b(x, A) = PQ.j{xEabX{x, A)), (138) 
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are the natural generalization of the 'squared solutions' introduced first for the 
s/(2)-case [32, 30]. By Pqj we have denoted the projector Pqj = ad J^ad j on the 
block-off-diagonal part of the corresponding matrix- valued function. 
The 'squared solutions' obviously satisfy the equation: 

i^ + [Q{x) - AJ, cabix, A)] = 0. (139) 

ax 

The indices a, b in (138) are taking values in appropriate ranges; for convenience 

below i, h, k (resp. by /, r, s) will be taking values in the ranges: 

1 < i,h,k < n, n+ 1 < l,r,s < n + m. (140) 

and by i < r we will mean 1 < i < n and n+l<r<n + m. We also introduce: 

where by eab we mean the relevant non- vanishing blocks of the matrices Eab', e.g.: 



^*"-(, OJ' ^'^ \Oers 

^ ^ ^ ' (142) 

F -{Oe^r\ p _{ 0' 

*"" ^0 J' "~ U" 0, 

Using these definitions and eqs. (28), (29) we find that the 'squared solutions' can 
be expressed in terms of the Jost solutions as follows: 

4(^, A) = a+(A)^+(x, A), e+(x,A) = c+(A)$+(x, A), 

4(^> A) = a+,(A)e+(x, A), e+(x,A) = c+ (A)H+(x, A). 

e,7(x,A) = aj;(A)^r^(x, A), 6,7 (x, A) = cr^(A)$^7,(x, A), 

e»~ft(a^' A) = c^kWQ^ki^, A), e" (x, A) = aj^{X)E;^{x, A). 

Here we assume summation over the repeated indices in the relevant range (140). 

We will also need the skew-scalar product: 

1 r°° 
[[X, Y]]=- dxtr {X{x), [J, Y{x)]) , (145) 

which is non-degenerate for X{x),Y{x) € M.. Then we can write down the 
reflection coefficients p^^^ in the form: 

P/+fe(A) = -i[[Q(y),e+(y,A)]]a+, (146) 

and similar expressions for the other reflection coefficients. Thus we have a for- 
mula analogous to the standard Fourier transform in which efi{y, A) can be viewed 
as the generalizations of the standard exponentials. 
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In order to work out the contributions from the discrete spectrum of L we will need 
the first two coefficients in the Taylor expansions of ^„5(x, A) and ^^i,{x, A); 



^a6(^> A) = ^%{x) + (A - \f)i>XAx) + 0{{\ - Xf)' 



^afe(^> A) = ^%{x) + (A - A±)6±^^.(x) + 0{{\ - Xff) 



where: 



$± . = $+ fx A±) $± . = "^^"^^""'^^ 



(147) 



(148) 



A=A- 



Indeed, from eq. (141) there follows that these are regular functions for all A G C±. 
However the 'squared solutions' e^j(x, A) have pole singularities of first order in 
the vicinities of A and therefore we will have 

4(^> A) =. ^^!^^ + 4^.(x) + 0{X - Xf), (149) 

where for a = i and b = I we have: 

and similar expressions for e^^^ (x) and e^^ (x) with different choices for the in- 
dices a, 6. 

The second type of Wronskian relations, which we will consider relate the variation 
of the potential SQ{x) to the corresponding variations of the scattering data. To this 
purpose we start with the identity: 



roo A 

lOO • ' , ^ 



X^x{x,X)\^^_^ = -i dx-—{ixSx){x,X). (151) 



oo 



dx 



To calculate the integrand in (151) we need to use the equation (133), satisfied by 
6x{x, A). Making use of eqs. (7) and (133) we get: 



r-OO 
lOO 



X'^x(a;,A)!;^_^ = i / dxx^Q{x)x{x,X). (152) 



oo 



We apply ideas similar to the ones above. Evaluating the l.h.side of (152) with 
x{x,X) ^X'^{x,X) andx(x,A) = x'{x,X) we find: 

v+^v+rr A)r -( «^'^«^(A) -ST+iX)c+iX)\ 

^ '^^ ^^'^^l«- l^<5/>+(A)a+(A) -c+5c+(A) j' ^^^^^ 

y-6y-(x Air = ( ^"^^"(^) 6p-{X)a-{X) \ 
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Multiplying by Eat and taking the trace we arrive at: 

<5p+(A) = 2i[[eUy,X),^d-/6Q{y)]]a^, 

5pT^{\) = -2i[[e+(y,A),ad7i<5Q(y)]]a+, 

These relations are basic in the analysis of the related NLEE and their Hamiltonian 
structures. Below we shall use them assuming that the variation of Q{x) is due to 
its time evolution, and consider variations of the type: 

5Q{x,t)=Qt5t + 0{{6tf). (156) 

Keeping only the first order terms with respect to 6t we find: 



dt 



2i[[eti{y^\)M'/Qi 



^^ = -2ih7(y,A),Q,(,,,.,,. 



(157) 



N 
' j-- 

N 



(159) 



(160) 



4.2. Completeness of the 'squared solutions' 

Let us introduce the sets of 'squared solutions' 

{*} = {*}cU{*}d, {*} = {*}cU{*}d, (158) 

{*}c = {*+(x,A), *r(x,A), i<r, A G m} , 

l^jd = "! ^riyV-^J) ^ri;j\-^)^ ^iryV"'"/' ir;j\''^) j ■_-, 

{*}c = {*+(x,A), *-,(x,A), i<r, Agr} 

i^/d = i ^jryA"'"/' ^ir;j\-^)j ^riyV-^J' ^ri;j\-^)j ._, 

where the subscripts 'c' and 'd' refer to the continuous and discrete spectrum of L. 
The 'squared solutions' in bold-face ^^■, ... are obtained from "^^■, ... by applying 
the projector Pqj, i.e. *+(x, A) = Poj^^^^^ix, A), see eq. (138). 

Theorem 3 (see [18]). The sets {^} and {^} form complete sets of functions in 
A4j. The corresponding completeness relation has the form: 

1 1"°° 
5{x-y)Iioj = - / dA(G+(x,y,A)-G-(x,y,A)) 

N (161) 

-2i^(G+(x,y)+G7(x,y)), 
where 

Hqj = Y,{Eir Eri - Eri Eir), (162) 
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G+{x,y,X) = ^e+(x,A)(8e+(y,A), 

i<r 
i<r 



(163) 



(164) 



t<r 



Idea of the proof: Apply the contour integration method to a conveniently chosen 
Green function, see [18]. D 

4.3. Expansions over the „squared" solutions 

Using the completeness relations one can expand any generic element F{x) of the 
phase space A4 over each of the complete sets of 'squared solutions'. We remind 
that F{x) is a generic element of A4 if it is a block-off-diagonal matrix-valued 
function, which falls off fast enough for |x| ^ oo. It can be written down in terms 
of its matrix elements F±{x) as: 

F{x) = Y,{F„{x)E„ + F„{x)Eri). (165) 



From (162) we get: 

-^tn ([J,F(x)] ® Il)noj = ^tr2noj(ll ® [J,F(x)]) = F{x). (166) 

where tr i (and tr 2) mean that we are taking the trace of the elements in the first 
(or the second) position of the tensor product. The result is: 

F{x) = - dXY, (ej>(a^, A)7^ .^(A) - e-{x,X)j^.^i{X) 



N 



(167) 



2i5](Z+^.(x) + Z^^^.(x)), 



i=i 



1 f°° / 

F{x) = — / dXY^ i(iri{x,X)-f^ i{X) - e-(x,A)7^i, 

TT J-00 ~Z. ^ 



N 



(168) 



where 

7|;a6(A) = [[e^a(y,A),F(y)]], (169) 
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Z+.{x) = Res Ve+(x,A)7^.,,(A), 
Zp.j{x) = Res_^e7.(x,A)7^.^.(A), 

'^'-'^'j i<r 

Z+.{x) = Res Ve+(x,A)7^.„(A), 
Z^.^.(x) = Res_ ^ e,r(a:, A)7^..,.(A), 



(170) 

_ _2^e„.(x,A)7p^^.(A), 

-^--^j i<r 



(171) 



'*'-'*'j i<r 

The completeness relation (161) is directly related to the spectral decompositions 
of A±, see e.g. [17, 25] for q = sl{n). It allows us to establish a one-to-one 
correspondence between the element F{x) G M and its expansion coefficients. 
Indeed, from eqs. (169) - (171) we can easily prove the following: 

Proposition 1, The function F{x) = if and only if all its expansion coefficients 
vanish, i.e.: 

7F;ir(^) = 7F;ri(^) = 0' ^ < ^5 Z+.ji^) = Z^.^.(x) = 0; (172a) 

7i^;„(A) = 7F;.r(A) = 0, i<r- Z+.^{x) = Z^.^.(x) = 0; (172b) 

where j = 1, . . . ,N. 

Proof. To show that from F{x) = there follows (172a) we insert F{x) = into 
the r.h. sides of the inversion formulae (169)-(171) getting (172). The fact that 
from (172a) there follows F{x) = is obtained by inserting (172a) into the r.h. 
side of (169)-(171). The equivalence of F{x) = to (172b) is proved analogously. 

4.4. Expansions of Q{x). 

Here we evaluate the expansion coefficients for F{x) = Q{x). As the reader have 
guessed already, their evaluation will be based on the Wronskian relations (136), 
(137) which we derived above. From them we have: 

[[e+(x,A),Q(x)]] = \dT^{\), [[e,:j(x,A),Q(x)]] = i4(A), (173) 

[[e+(x,A),Q(x)]]=ib+(A), [[e~{x,X),Q{x)\]=\bT^{\), (174) 

Skipping the calculational details we get the following expansion of Q{x) over the 
systems {$^} and {*^}: 

Q{x) = - dA5](r+(A)*+(x,A)-T,7(A)*7,(x,A 

N (175) 

k=l i<r 
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N (176) 

fc=l i<r 

4.5. Expansions of ad ^j^5Q{x). 

Here we evaluate the expansion coefficients for F{x) = ad 'j^5Q{x). Their evalu- 
ation is based on the Wronskian relations (153), (154). Now we have: 



[[e+(x,A),ad7i5Q(x)]] = '-{{5t+)c+ {\)\u 
[[e-^{x,\)M-/mx)]] = -\{{6t-)c-{\))u. 



B,7(x,A),ad7i5Q(x)]] = -{{5p-)a-{\))a, 



(177) 



(178) 



Skipping the calculational details we get the following expansion of ad j 6Q{x 
over the systems {$^} and {*^}: 

^AfbQix) = — y_ dAE (5t+(A)*+(x,A) + 5t,-t(A)*;,(x,A 

°° i<r 

N (179) 

+ EE('^'^r;,(^)-^X7;,W)' 
fc=l i<r 

adfSQix) = ^|_^dAE('5/>i(A)*+(x,A) + 5p-(A)*-(x,A 

^ '^' (180) 

k=l i<r 

where 

^'wX,i^) = 6XfTX,^%-Ax) + 'J4^*i,(^), (181) 

S'wX,ix) = 6XfpX,^X^{x) + 5p± ^,*^,;,(x) (182) 

The expansions (175), (176) combined with proposition 1 is another way to estab- 
lish the one-to-one correspondence between Q{x) and each of the minimal sets of 
scattering data 7i and T2 (56). Likewise the expansions (179), (180) and proposi- 
tion 1 establish the one-to-one correspondence between the variation of the poten- 
tial 6Q{x) and the variations of the scattering data 6T1 and 6T2. 



BASIC ASPECTS OF SOLITON THEORY 33 

4,6. The generating operators 

To complete the analogy between the standard Fourier transform and the expan- 
sions over the 'squared solutions' we need the analogs of the operator Dq = 
— id/dx. The operator Dq is the one for which e'^^ is an eigenf unction: Dqg^^^ = 
^giAx xherefore it is natural to introduce the generating operators A± through: 

(A+ - A)*+(x,A) = 0, (A+ - A)*-(x, A) = 0, 

(A_ - A)*+(a;, A) = 0, (A_ - A)*;^(x, A) = 0. 

Their derivation starts by introducing the splitting: 

ea6(^>A) = e^f (x,A) + e±(x,A), e^f (x,A) = (11 -Poj)e^,(x, A), (184) 

into the equation (139). Then eq. (139) splits into: 

i-^ + [Q(x),e±(x,A)] =0, (185) 

i^ + [Q(.x),etl^{x,X)] = A[J,e±(x,A)], (186) 

Eq. (185) can be integrated formally with the result 

e^f (x, A) = ctl%W+\r dy [g(y),e± (y, A)], (187) 

Next insert (187) into (186) and act on both sides by ad j^. This gives us: 

(A±-A)e±(x,A) = i[C,^^^(A),ad7iQ(x)], (189) 

where the generating operators A± are given by: 

(190) 



A±X(x) = a.dj^[ i^ + i 



dX 
dx 



Q(x), / dy[Q{y),X{y)] 

itoo 



Thus e^j(x, A) will be an eigenfunction of A+ or A_ if only if C^^ {y, A) = 0. 
Evaluating the limits of (188) for all combinations of indices a,bv/e find (i < r): 

(A+ - A)*+(x, A) = 0, (A+ - A)*r (x, A) = 0, 

(A+ - A+)*,+ ^.(x) = 0, (A+ - AT)*r ^.(x) = 0, 

(A. - A)*+(x, A) = 0, (A_ - A)*;,(x, A) = 0, 

(A_ - A+)*+^^.(x) = 0, (A_ - AT)*- ^.(x) = 0, 

The rest of the squared solutions are not eigenfunctions of neither A+ nor A_ : 

(A+ - X+)4fti.j{x) = *+.^-(x), (A+ - A7)*,;./x) = ^;,.j{x), 

(A_ - A+)*t;i(a^) = *»>;j(2;), (A- - A-)*riy(x) = ^-^(x), 
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i.e., ^^j (x) and ^^^ (x) are adjoint eigenfunctions of A_|_ and A_. This means 
that Xf,j = 1, . . . , A^ are also the discrete eigenvalues of A± but the correspond- 
ing eigenspaces of A± have double the dimensions of the ones of L; now they 

are spanned by both ^ab-ji^) ^^^ ^a6;i(^)- Thus the sets {^} and {<!>} are the 
complete sets of eigen- and adjoint functions of A^, and A_ . 



5. Fundamental Properties of the MNLS-type equations 

In this Section we describe the fundamental properties of the NLEE. 

5.1. The class of the MNLS-type equations 

Let us insert variations of the form (156) into the expansions (179), (180). This 
results in: 



idQ i r ^xV-f dT+ , , ^, , dr, 



ad 

-CXD 



■^ dt 2ttJ- 



i/>^i:Hf*J(-^)+^*'<(-^> 



i<r 

N 

fc=l i<r 



(194) 



i<r \ 

N 



(195) 



where 



fc=l i<r 



.',±.. .^_d^^.± ^± r^^ , ^^i^± 



W::,%{x,t) = ^r,t;^.*^,^^.(x) + ^$±^^.(x), (196) 

Next from (191), (192) there follows that: 

(/o(A+) - /o(A))*+(x, A) = 0, (/o(A+) - /o(A))*- (x, A) = 0, 

(198) 
(/o(A+) - /o(A+))*+ ^.(x) = 0, (/o(A+) - /o(A7))*,r;/x) = 0, 

(/o(A+) - /o(A+))*,+,,.(x) = /(A+)*+ ^.(x), 

(/o(A+) - /o(A-))*i,y(x) = t\Xj)^T,,j{x), 
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and similar relations between A_ and its eigenf unctions. Combining them with the 
expansions (175), (176) we get: 

/o(A_)Q(x) = - / dA/o(A)^ (t+(A)*+(x,A) -T,7(A)*;,(x,A) 






N (200) 

+ ^T.I1 (/o(A+)r+,*+;,(x) + uxjK'i.,,^;r,ji^)) , 

k=l i<r 

/o(A+)Q(x) = — / dXj2foiX) (/>+(A)*+(x,A) -p,r(A)*,r(x,A) 



(201) 



2 E E (/o(A+)p+y*+;,(x) + /o(A^. )Pir.j'^„-A^: 



fc=l i<r 

Now we can prove that the principal series of MNLS-type equations has the form: 

mdf^ + fo{A)Q{x,t) = 0, (202) 

where A can be either A_|_ or A_ and /o(A) determines the dispersion law of the 
corresponding NLEE. 



Theorem 4. The NLEE (202) are equivalent to each of the following evolution 
equations for the scattering data of L: 

i^T2/o(A)p± = 0, ^ = 0, i^^2/o(A±),J = 0, (203) 
i^±2/o(A)r± = 0, ^ = 0, .^±2/o(A±)rJ = 0, (204) 

Proof: Insert the expansions (195) and (201) into the left hand side of the NLEE 
(202) and use proposition 1. This immediately proves the equivalence of the NLEE 
to the linear equations (203). Analogously, from the expansions (194), (200) and 
proposition 1 one proves the equivalence between the NLEE and (204). D 

5.2. Integrals of motion - principal series. 

Dealing with the block Zakharov-Shabat system we can distinguish two types of 
NLEE and, consequently, two types of series of conservation laws that provide the 
Hamiltonians of the NLEE. 

The MNLS equations have maximally degenerated dispersion law and belong to 
the principle series of NLEE. Their Hamiltonians have local densities, i.e., their 
densities depend on Q{x,t) and its x-derivatives. At the same time they have 
maximal number of generating functionals of conservation laws: the whole blocks 
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a^(A) and c^(A). However not all of these functional are in involution, see 
Subsection 5.4 below. 

The NLEE characterized by generic dispersion laws have as a rule non-local Hamil- 
tonian densities and a minimal possible number of generating functional of inte- 
grals of motion. These are provided by the invariants (the eigenvalues) of a^(A) 
and c^(A). Special combinations of these as in eq. (61) produce the Hamiltonians 
of the MNLS type equations. In this Subsection we will concentrate on those. 
In Subsection 2.6 we showed how these functionals A'^(A) and C^(A) can be 
expressed through the minimal sets of scattering data, see eq. (62). Our aim here is 
to derive a recurrent procedure which would allow one to express their coefficients 
in the asymptotic expansions: 

OO CO 

A+{X) = J2 h\'\ C+(A) = Y. Ju\~\ (205) 

k=l k=l 

in terms of (5(x, t). To this end we make use of a Wronskian type relation: 

(11 - Poj) x^{x, A) (x+ + ix Jx+) ir=_oo = ( "^o"^ -£+c+ ) ' (2^^^ 

valid for A G C+. Here and below by 'dot' we denote derivative with respect to A. 
In what follows we will need also the standard formula: 

tr In a+ ( A) = In det a+ ( A) , (207) 

and its consequence 

— tr lna+(A) = In det a+(A)a+(A). (208) 

dA 

Next we express x^i^^ ^) ^^ the form: 

X+{x, A) = (11 + W{x, X))Z{x, A)e-^^-'^, (209) 

where 

W(xX)-( ° W('\x,X)\ zfxA)-(^^^'^(^'^) ° ^ 

(210) 
Inserting (209) into (7) and separating the block-diagonal and block-off-diagonal 
parts we get for W{x, A) and Z{x, A) the following system: 

dW 

i^ + Q{x,t)-WQW{x,t,X) = X[J,W{x,t,X)], (211) 

ax 

A 7 

i--Z(x, t, A) + Q(x, t)W(x, t, A) = 0. (212) 

ax 
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Combining eqs. (209) and (206) we obtain: 

_ /a+(A)d+(A) 



Z{x,t, X)Z 







-c+(A)c+(A)]' 



(213) 



which leads, in view of (208) and (212) to 



Z{x,t, X)Z 



oo 

X= — OD 



i dxZ{x,t,X)Q{x,t)W{x,t,X)Z{x,t,X). (214) 



If we multiply both sides of (214) by ( p, p, I and ( p. ,, 1 take the trace and inte- 



11 



II 



(215) 



grate over A we get: 

yl+(A) =lndeta+(A) = i H dxtr (Q{x,t)W{x,t,X) ( jj 

C+(A) =lndetc+(A) = -i H dxtr (Q{x,t)W{x,t,X) (^ 
Using the asymptotic expansions of yi+(A), C+(A) (see eq. (205)) and W{x, t, A): 

oo 

Wix,t,X) = J2Wkix,t)X-\ (216) 

fe=i 

we arrive at the following expressions for the integrals of motion Ik and J^ : 

II 



11 



Jk 



— oo 
oo 



dxtr {Q{x,t)Wk{x,t) 
dxtr (Q{x,t)Wk{x,t) 



(217) 



The last step in these considerations consists in deriving recurrent relations for 
calculating Wk{x,t) in terms of Q{x,t) and its x-derivatives. This is done by 
inserting (216) into the equation (211) with the result: 



W,{x,t) = ^[J,Q{x,t)], 



Wk+iix,t) = - 



In particular we have: 



J,i^- E Wp{x,t)Q{x,t)Ws{x,t) 

p+s=k 



dx 



(218) 



W2{X, t) = -Q^, W3(X, t) = — \j, Q^^ + Q3 

4 lb L 



(219) 



which leads to 

^oo 



J ("OO ^ re 

h = - dxtr(gr), Ji = -- / 

'^ J — OO 41: J — ( 



dxtr (rg), 
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1 /-OO 1 /-oo 

h = --7 dxtr (qr^), J2 = - / dxtr (rq^), (220) 

4 J— 00 4 J —00 

[ roo i /■°° 

/a = -- / dz tr (gr^^a; + gr^r), -^3 = -q / dxtr (rg^^^a; + rgrg'). 

O J-00 o J-00 

5.3. Hamiltonian properties of the MNLS eqs. 

Let us briefly outline the Hamiltonian properties of the NLEE (202). Obviously the 
MNLS describes an infinite dimensional Hamiltonian system with Hamiltonian: 

i^^MNLS = 7: dytr [QQ^^ + Q^ix,t)) = 4i{h + ^3), (221) 

and Poisson brackets: 

{QksiVi t),rri{x, t)} = iSik5rsS{x - y), (222) 

or, equivalently, by the canonical symplectic form: 

/•OO 

r^o = i / dxtr {5r{x) A5q{x)) 



00 



1 /""^ / \ 

= - / dxtr (ad j^6Q{x) A [J, ad 'j^6Q{x)) . (223) 

The second expression is preferable to us because it makes obvious the interpreta- 
tion of SQ{x, t) as local coordinate on the co-adjoint orbit passing through J. It is 
also expressed through the skew-scalar product by: 

VLq = \ [[ad -/5Q /\ ad-/ 5Q^. (224) 

It can be evaluated in terms of the scattering data variations. To do this we insert 
the expansion of ad~/5Q into VLq and then use again the equations (177), (178). 
After some calculations we get: 

^0 = ^ j^ d\ (0^(A) - n,{\)) - i Y. {^l^ + f^o;^.) , (225) 



-ti ^^u/ I. ,^u^J^ yy\)) , ^lqj 



J]±(A) = -tr {6t'^c'^ a (5/9^a^(A)) , Q^- = Res n^{X), (226) 



Here we skip the explicit expressions for JIq which are rather involved. From 
(225) it is not even obvious that JIq is closed. 

The Hamiltonian formulation of the MNLS eq. with Qq and Hq is just one member 
of the hierarchy of Hamiltonian formulations of MNLS provided by: 

nk = j[[adf 6Q A A'' ad f5Q]], A=^(A+ + A_), (227) 
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Hk = 4i(4+3 + Jk+s). (228) 

We can also calculate Q^ in terms of the scattering data variations. Doing this we 
will need also eqs. (198), (199). The answer is 

1 /'CXD 

^k = ^J dxx''{n+{x)-n,{\))-iJ2{n+^ + n~.^^), (229) 

nti = Res X''Q^{X). (230) 

This allows one to prove that if we are able to cast Qq in canonical form then all 
r^fc will also be cast in canonical form and will be pair-wise equivalent. 

5.4. The classical i?-matrix and the NLEE of MNLS type 

One of the definitions of the classical i?-matrix is based on the Lax representation 
for the corresponding NLEE. We will start from this definition, but before to state 
it will introduce the following notation: 

^Uix,X)(^U{y,ii)Y (231) 

which is an abbreviated record for the Poisson bracket between all matrix elements 

of U{x, X) and U{y, fi) 

\uix,X)(i^U{y,ii)\ ={Uik{x,X),Uimiy,f^)}. (232) 

L ) ik,lm 

In particular, iiU{x, A) is of the form: 

U{x, A) = Q{x, t) - A J, Q{x, t) = ^{qirEir + PriEri), (233) 

i<r 

and the matrix elements of Q{x, t) satisfy (222), then: 

\u{x,X)®U{y,^i)\ = iJ2(^ir Eri - Eri ® Eir)5{x - y). (234) 

^ ^ i<r 

The classical it!-matrix can be defined through the relation [12]: 



U{x, A) U{y, /u) ^ = i [R{X - n), U{x, A) (g) II + II ® U{y, ^)] 6{x - y). 

(235) 
which can be understood as a system of N'^ equation for the N'^ matrix elements 
of R{X — /[/). However, these relations must hold identically with respect to A and 
H, i.e., (235) is an overdetermined system of algebraic equations for the matrix 
elements of R. It is far from obvious whether such R{X — fi) exists, still less 
obvious is that it depends only on the difference X — ii. In other words far from any 
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choice for U{x, A) and for the Poisson brackets between its matrix elements allow 
i?-matrix description. Our system (235) allows an i?-matrix given by: 

^(A-m) = -^Y^' ^236) 

where P is a constant N"^ x N"^ matrix: 

TV 

P= J2 ^-b^ ^ba- (237) 

a,b=l 

The matrix P possesses the following special properties: 

P{X®Y) = {Y(S)X)P, P^ = i, (238) 

i.e., it interchanges the positions of the elements in the direct tensor product. By 
using these properties of P we are getting: 

[P, Q{x) «) II + II ^ Q{x)] = 0, (239) 

i.e., the r.h.side of (235) does not contain Q{x, t). Besides: 

[P, A J ® II + ^11 (g) J] = (A - ^) [P, J (g) II] 

= -2(A - li) \y^{Eir ® Eri - Eri ® Eir) J . (240) 

\i<r / 

where we used the commutation relations between the matrices Eah'- 

[Eab, Ecd] = Eadhc — Ech^da- (241) 

The comparison between (239), (240) and (235) leads us to the result, that R{X—fi) 
(236) indeed satisfies the definition (235). 

Let us now show, that the classical P-matrix is a very effective tool for calculating 
the Poisson brackets between the matrix elements of r(A). It will be more con- 
venient here to consider periodic boundary conditions on the interval [—L, L], i.e. 
Q{x — L) = Q{x + L) and to use the fundamental solution T{x, y, A) defined by: 

.dT{x,y,X) 



dx 



+ Uix,\)T{x,y,X) =0, T{x,x,X) = i. (242) 



Skipping the details we just formulate the following relation for the Poisson brack- 
ets between the matrix elements of T{x, y, A): 



T{x, y, A) T{x, y,fJ.)y = [P(A - fi),T{x, y, A) O T{x, y, /i)] 



(243) 



The corresponding monodromy matrix Tl{X) describes the transition from — P to 
L and Tl{X) = T{—L, L, A). The Poisson brackets between the matrix elements 
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of Tl{\) follow directly from eq. (243) and are given by: 

{tl{\)®Tl{A = [R{X-fi),TL{X)(^TL{fi)]. (244) 

An elementary consequence of this result is the involutivity of the integrals of 
motion /^ ^ from the principal series which are from the expansions of: 

CO CO 

Indet a+(A) = J2 lL,kX~\ - Indet c^(A) = ^ lL,k\-\ (245) 

k=l k=l 



oo oo 



Indet c+(A) = J2 JL,k\-\ -Indet a^(A) = ^ JL,k)^^\ (246) 

fc=i fc=i 

An important property of the integrals I^^ and J^^ is their locality, i.e. their 
densities depend only on Q and its z-derivatives. 

The simplest consequence of the relation (243) is the involutivity of lL,k, JL,k- 
Indeed, taking the trace of both sides of (243) shows that {tr Ti(A), tr Ti^fi)} = 0. 
We can also multiply both sides of (243) by C ® C and then take the trace using 
eq. (239); this proves: 

{trTL(A)C,trrL(^)C} = 0. (247) 

In particular, for C = II + J and C = II — J we get the involutivity of: 

|tra+(A),tra+(^)| = 0, ftra^(A), tr a^(^)| = 0, 

|trc+(A),trc+(^)| = 0, |trc^(A), tr c^(^)| = 0, 

Eq. (243) was derived for the typical representation V^'^^ of (S ~ SU{n + m), 
but it holds true also for any other finite-dimensional representation of (S. Let 
us denote by V'^^'* ~ a'^V'^^^ the /c-th fundamental representation of ©; then the 
element Tl{\) will be represented in V^^> by a'^Tl(A) - the k-th. wedge power of 
Tl{\), see [28]. In particular, if we consider eq. (243) in the representation y(") 
and sandwich it between the highest and lowest weight vectors in !/'-"•* we get: 

{deta+(A),deta+(^)} = 0, {det c^(A),det c^(^)} = 0. (249) 

Likewise considering (243) in the representation y(™) and sandwich it between 
the highest and lowest weight vectors in y('") we get: 

{deta^(A),deta^(^)} = 0, {det c^(A),det c^(^)} = 0. (250) 

Since eqs. (249) and (250) hold true for all values of A and ^ we can insert into 
them the expansions (245) with the result: 

{/L,fc,/L,p} = 0, {JL,fc,JL,p} = 0, A:,p = l,2,.... (251) 
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Somewhat more general analysis along this lines allows one to see that only the 
eigenvalues of a^(A) and c^(A) produce integrals of motion in involution. 

Taking the limit L ^ cx) we are able to transfer these results also for the case 
of potentials with zero boundary conditions. Indeed, let us multiply (243) by 
E{y, A) E{y, /x) on the right and by E^^{x, A) E^^{x, fi) on the left, where 
E{x, A) = exp(— lA Jx) and take the Hmit for re ^ oo, y ^ — oo. Since: 

gia;{A-/i) 

lim — = ±[tt5{X - fi), (252) 

x^iboo A — H 

we get: 

TiX)0T{fi)\ = R+{X - fi)T{X) <g) Tip) - T(A) (g) r(^)ii_(A - fi), 
R±iX - ij) (253) 

= -TwT r \^Eik(^ Eki +^Ers^ Esr ± ivr(5(A - n)Iloj, 

^^^ ~^^) \ik rs J 

where IIoj is defined by eq. (162). Analogously we prove that: 

i) the integrals I^ = lim/,^oo lL,k and Jp = YvaiL^ao Jl,p are in involution, i.e.: 

{h, Ip} = {h, Jp} = {Jk, Jp} = 0, 

for all positive values of k and p; ii) only the eigenvalues of a^(A) and c^(A) 
produce integrals of motion in involution. 

5.5. Generic NLEE 

Now we apply the expansions above to the analysis of the generic NLEE related to 
L. Each of these NLEE is determined by its dispersion law: 

The analysis of these NLEE is based on the expansions [Bk, ad j^Q{x, t)] over the 
systems {*} and {*}. The corresponding expansion coefficients are obtained by 
multiplying eq. (136) on the right by B^ and taking the trace. Then apply A*^ to 
both sides of the expansions. This proves a theorem generalizing theorem 4. 

Theorem 5. The generic NLEE with polynomial dispersion law /(A) (254) are of 
the form: 

iadj^-^ + J2^''[Bk,ad-]^Q{x,t)] =0, (255) 

•^^ fc=i 



BASIC ASPECTS OF SOLITON THEORY 43 

and are equivalent to each of the following evolution equations for the scattering 
data of L: 

i^ - /2(A)p+(A, t) - p+(A, t)/i(A) = 0, 



dt 
dp' 



■■"^^ +/i(A)p-(A,t)+p-(A,t)/2(A) = 0, 



(256) 



± 






+ /up7 it) - p7 {t)f2,j =0, ^ = 0' 

+ /i(A)r+(A,t) + r+(A,t)/2(A) = 0, 
-/2(A)r-(A,t)-r-(A,t)/i(A) = 0, 



dt 
dr+ 
i-TT + ftA^) + ^/(*)/2!'j = 0, /!%• = /l,2(A,± 



(257) 



dt 



drr dAt 

' -^ - fi,r- (t) - T- it)fr^ = 0, -rf = 0. 



dt ■'^'■^ -^ ' ' J ' '■'^^J ' dt 

As a consequence of theorem 5 we get that S^{X) and T^(A) satisfy Unear set of 
equations: 

d <?=•= dT^ 

i^ + [/(A),5±(A,t)] =0, i^ + [/(A),T±(A,t)] = 0, (258) 

which can be easily solved. In addition the block-diagonal matrices D^{X) are not 
integrals of motion but rather satisfy: 

d D^ 
i^^+[/(A),D^(A,t)] =0, (259) 

From eq. (259) there follows that for generic /(A) only the invariants of D^{X, t) 
(or, equivalently, only the invariants of a^(A, t) and c^(A, t)) provide series of in- 
tegrals of motion in involution. As a consequence the generic NLEE (255) possess 
soliton solutions whose velocities may depend on time. Examples of such NLEE 
and the properties of their soliton solutions called boomerons and trappons have 
been analyzed by Calogero and Degasperis [7, 8, 9, 10]. 

6. Discussion 

In order to understand better the idea that the expansions over the 'squared solu- 
tions' are generalized Fourier transforms we will outline the limit to small poten- 
tials, i.e. Q{x) ^ 1, for which the Born approximation is adequate. We will see 
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that in this Umit all these expansions turn out to be just the usual Fourier trans- 
forms. 

It is known [1] that the corresponding Zakharov-Shabat system may have discrete 
eigenvalues only provided Q{x) is large enough, while in the Born approxima- 
tion it does not possess discrete spectrum. The eigenfunctions of the continuous 
spectrum are well approximated by the 'plane waves': 

X+(x,A)~x"(x,A)~e-i^-^-, (260a) 



and the scattering data are provided by the Born approximation: 

)e-2i^^ dx 

(260b) 



fOO 

pIW = i / g^r(a;)e''^" dx, p;,{\) = i / rri{x)e'^''^- dx. 



oo 



a+(A) ~ 1, a-(A)~l, 

As a result the "squared solutions" are also approximated by 'plane waves': 

*+(x. A) ~ *7,(x, A) ~ Eire^''^^ *r (x, A) ~ *+(x. A) ~ E^e-^''^^ 
The completeness relation (161) acquires the form: 

5{x - y)noj = - / dx V (e-2'^("-^)s,, Eri - e^'^^""^)^,, E, 



of the usual Fourier transform for matrix-valued functions. The recursion operators 
K± in this limit go into the purely differentiation operator: 

d 



Af • ~ - 
* 4 



' dx 



(261) 



whose spectrum is purely continuous. 

On the other hand the expansions over the "squared solutions" can be understood as 
expansions over the eigenfunctions of A+ and A_. Note, that while id/dx does not 
have discrete spectrum, the operators A± may have discrete eigenvalues whenever 
L has. Finally, if we restrict ourselves to the special class of reflectionless poten- 
tials for which p^(A) = for all A G M, we obtain finite dimensional subspace of 
M. Then the operators A± operate nontrivially only in this subspace. 

7. Conclusions 

We showed that the interpretation of the ISM as a generalized Fourier transform 
holds true for the generalized Zakharov-Shabat systems related to the symmetric 
spaces SU{n + m)/ S{U{n) ® U{n)). Using the standard dressing method we 
outlined the structure of the singularities of the fundamental analytic solutions and 
constructed the soliton solutions with multidimensional eigenspaces (the projectors 
Pj (x) may have rank Vj > 1). 
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The expansions over the 'squared solutions' are natural tool to derive the funda- 
mental properties not only of the MNLS type equations, but also of the NLEE with 
generic dispersion laws. Some of these equations, besides the intriguing proper- 
ties as dynamical systems allowing for boomerons, trappons etc., may also have 
interesting physical applications. 

Another interesting area for further investigations is to study and classify the re- 
ductions of these NLEE. For results along this line for the MNLS equations see the 
reports [21] and [22]; reductions of other types of NLEE have been considered in 
[16, 19, 20, 33]. 

One can also treat generalized Zakharov-Shabat systems related to other symmet- 
ric spaces. This would require substantial changes in the dressing factors. The 
expansions over the 'squared solutions' can be closely related to the graded Lie 
algebras, and to the reduction group and provide an effective tool to derive and 
analyze new soli ton equations. For more details and further reading see [11, 16]. 
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